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We calculate the nuclear matrix elements of the neutrinoless double beta (Oi^/3/3) decays of ^^Ge 
and ^■^Se for the light-neutrino exchange mechanism. The nuclear wave functions are obtained by 
using realistic two-body forces within the proton-neutron quasiparticle random-phase approximation 
(pnQRPA). We include the effects that come from the finite size of a nucleon, from the higher-order 
terms of nucleonic weak currents, and from the nucleon-nucleon short-range correlations. Most 
importantly, we improve on the presently available calculations by replacing the rudimentary Jastrow 
short-range correlations by the more advanced unitary correlation operator method (UCOM). The 
UCOM corrected matrix elements turn out to be notably larger in magnitude than the Jastrow 
corrected ones. This has drastic consequences for the detectability of QvPfi decay in the present and 
future double beta experiments. 
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The nuclear matrix elements of neutrinoless double 
beta (0v(3(3) decay have become an important issue in the 
present-day neutrino physics (see, e.g., P, 3, This 
has been boosted by the verification of the existence of 
neutrino mass by the oscillation experiments Q and the 
claimed discovery of the Ovpf) decay 0|. Further in- 
centive to produce reliable nuclear matrix elements comes 
from the needs of the running NEMO 3 [i| and CUORI- 
CINO experiments, as well as the future large-scale 
experiments under R&D planning and construction (see, 
e.g., For these experiments the nuclear matrix el- 

ements are an essential prerequisite in extraction of re- 
liable values for the absolute mass scale of the neutrino 
fill and possibly the CP phases of the neutrino-mixing 
matrix 

The most popular nuclear model to treat the struc- 
ture of medium-heavy and heavy double-beta decaying 
nuclei is the proton-neutron q uasi par ticle random-phase 
approximation (pnQRPA) [id. Il4 |l 5|. Also some recent 
shell-model results are available [16|, 1131 . The pnQRPA 
is a model tailored to describe, in an efficient way, the 
energy levels of odd-odd nuclei and their beta decays to 
the neighboring even-even nuclei [l3|. Also its deriva- 
tive, renormalized pnQRPA [l^, has been used (20l. [2l| 
to compute double-beta matrix elements, although its use 
has been heavily criticized (see, e.g., [22| and references 
therein). A problem with the use of the pnQRPA (and 
the renormaHzed pnQRPA) is that it contains a free pa- 
rameter, the so-called particle-particle strength param- 
eter, jpp, that controls the magnitude of the proton- 
neutron two-body interaction matrix elements in the 
T = pairing channel [13, [131 ■ There are basically two 
ways to fix the value of this parameter, either by using 
the data on two-neutrino double beta. (2iy00) decay [2l| 
or the data on single beta decay [25l. l2a|. In the case of 
^^Ge and ^^Se there is no available data on single beta 
decays so that in this work we have chosen to use the 
2j//3/3 data to fix the possible values of gpp. 

In this article we address the mass mode of the Oi^/3/3 



decay where a light virtual Majorana neutrino is ex- 
changed by the two decaying neutrons of the initial nu- 
cleus. Typically the exchanged momentum is so large 
as to force the two neutrons to overlap unless steps are 
taken to prevent the occurance of such a spurious event. 
The traditional way (27l . to remove this spuriosity 
is to introduce an explicit Jastrow type of correlation 
function into the involved two-body transition matrix el- 
ements in the parametrization of Miller and Spencer [29| . 
This method, although microscopically inspired, is just a 
phenomenological way to introduce short-range correla- 
tions into the two-nucleon relative wave function. A con- 
spicuous fiaw of the Jastrow method is that the Jastrow 
function effectively cuts out the small r |3fl|| part from 
the relative wave function of the two nucleons. For this 
reason, the traditionally adopted Jastrow procedure 
does not conserve the norm of the relative wave function 

m. 

In the present calculations we improve on the Jas- 
trow method and adopt the more sophisticated micro- 
scopic approach of unitary correlation operator method 
(UCOM) Q. In the UCOM one obtains the correlated 
many-particle state \^) from the uncorrelated one as 



(1) 



where C is the unitary correlation operator [32|. Due 
to the unitarity of the operator C, the norm of the cor- 
related state is conserved and no amplitude is lost in 
the relative wave function. In the calculations this 

leads to a more complete description of the relative wave 
function at small distances r, as was demonstrated in 
[33 |. It should be stressed that no extra free parame- 
ters are introduced by the use of the UCOM at the level 
of double-beta-decay calculations. All the needed pa- 
rameters have been fixed by minimization procedures for 
/3;9-independent observables [34]. The UCOM method 
has been demonstrated [35j to produce good results for 
the binding energies of nuclei over a wide mass range al- 
ready at the Hartree-Fock level. It was also shown that 
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the UCOM renders a good starting point for inclusion of 
long-range correlations by means of many-body pertur- 
bation theory. 

In [s^l it was demonstrated for the ^^Ca and ""^Ge 0^(3(3 
decays that the Jastrow procedure leads to the excessive 
reduction of 30% - 40% in the magnitudes of the Oi^/3/3 
nuclear matrix elements. At the same time the UCOM 
reduces the magnitudes of the matrix elements only by 
7% - 16%. This explains the large short-range correlation 
corrections to the matrix elements of (2l| . 

The issue about the magnitude of the short-range cor- 
rections is an extremely important one since it directly 
affects the magnitudes of the relevant nuclear matrix el- 
ements used to extract the neutrino masses from poten- 
tially succesful future double beta experiments. There 
are large differences between the Jastrow and UCOM 
corrections, the Jastrow corrected matrix elements be- 
ing substantially smaller that the UCOM corrected ones. 
This difference would severely alter the predicted sensi- 
tivities of future neutrino experiments, the UCOM cor- 
rected matrix elements being more favourable for the de- 
tection of neutrinoless double beta decay. 

The double beta decays of ''^Ge and ^^Se proceed 
through the virtual states of the intermediate nuclei ""^ As 
and *^Br to the ground states of the final nuclei ^^Se and 
^^Kr. By assuming the neutrino-mass mechanism to be 
the dominant one, we can write the inverse of the half-life 
as fil 
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where rric is the electron mass and cf^^ is the leptonic 
phase-space factor. The Oi//3/3 nuclear matrix element 
TVf (O'') consists of the Gamow-Teller, Fermi and tensor 
parts as 



(2) 



with a standard parametrization, optimized for nuclei 
near the line of beta stability. Slight adjustments were 
done for some of the energies at the vicinity of the pro- 
ton and neutron Fermi surfaces to reproduce better the 
low-energy spectra of the neighboring odd-A nuclei and 
those of the intermediate nuclei. 

The Bonn-A G-matrix was used as a two-body inter- 
action and it was renormalized in the standard way, as 
discussed e.g. in Refs. [stI. fs^. Due to this phenomeno- 
logical renormalization we did not perform an additional 
UCOM renormalization [34] of the two-body interaction. 
After fixing all the Hamiltonian parameters the only free 
parameter left was the ^pp parameter mentioned earlier. 
In Fig. [T] we have studied the gpp dependence of the ma- 
trix element M(o^) of Eq. ^ for both ™Ge and ^^Se. 
We have used R = 1.2A^/^ im as the nuclear radius, and 
the finite-size, higher-order-term and UCOM corrections 
were taken into account. Here one can see the typical 
break-down of the pnQRPA at large values of gpp. More- 
over, the corresponding break-down points for the two 
nuclei are close to each other. 
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Numerical calculations show that the tensor part in 
([3]) is quite small and its contribution can be safely ne- 
glected in what follows. The expressions for the phase- 
space factor, the effective neutrino mass {ruv) and the 
matrix elements of ^ are given, e.g., in [ill, 113, [T3 |. 
To the "bare" matrix elements we have applied the Jas- 
trow short-range correlation corrections, together with 
the higher-order terms of nucleonic weak currents and 
the nucleon's finite-size corrections in the way described 
in [2ll. |36|. In addition, we have computed the corrected 
matrix elements by replacing the Jastrow correlations by 
the UCOM correlations. 

We calculated the wave functions of all the nuclear 
states in the intermediate nuclei by the use of the pn- 
QRPA framework in the model space of lp-0f-2s-ld-0g- 
Ohii/2 single-particle orbitals, both for protons and neu- 
trons. The single-particle energies were obtained from 
a spherical Coulomb-corrected Woods-Saxon potential 



FIG. 1: The OvPP nuclear matrix elements M^""^ of Eq. Q 
for the decays of ^^Ge and *^Se as functions of the particle- 
particle interaction parameter gpp. 



We obtained the physical values of gpp by using the 
method of [2l|. Consequently, we used the recommended 
data (lil on 2i^/3/3-decay half-lives of ^^Ge and ^^Se by 
including the experimental error limits and the uncer- 
tainty in the value of the axial-vector coupling constant 
1.0 < gA < 1.254. The resulting intervals [40|| of "ex- 
perimental matrix elements" were then converted to the 
following intervals of gpp values 



1-02 < 5pp < 1.06 for ^^Ge 
0.96 < 5pp < 1.00 for ®^Se . 



(4) 



In Fig.[2]we show for the *^Se decay the decomposition 
of the total matrix element ^ into multipoles. The bare 
matrix element contains no short-range, finite-size or 
higher-order-term corrections, whereas the UCOM and 
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FIG. 2: Multipole decomposition of the matrix element M^""' 
of Eq. ^ for the decay of *^Se. The spread in the multipole 
contributions corresponds to the Qpp interval for *^Se in ([4]). 



TABLE I: Matrix element M^""' of Eq. ^ computed by cor- 
recting successively the bare matrix element (b.m.e.) by the 
higher-order terms in the nucleonic current (A), by the nu- 
cleon finite-size effect (B), and by either the Jastrow (C) or 
UCOM (D) correlations. The values ^pp = 1.0 and = 1-254 
were used in the calculations. 



Nucleus 


b.m.e. 


+A 


+A+B 


+A+B+C 


+A+B+D 


'"^Ge 


-8.529 


-7.720 


-6.356 


-4.723 


-6.080 


«^Se 


-5.398 


-4.862 


-3.914 


-2.771 


-3.722 



Jastrow matrix elements include all these corrections. 
The spread in the multipole contributions corresponds to 
the gpp interval for ®^Se in More specifically, the up- 
per end of the bar represents the case gpp — 0.96) with 
gA = 1-0 and the lower end the case gpp — 1.00 with 
gA = 1.254. In the figure we see that the corrections 
substantially reduce the magnitude of a given multipole 
contribution. The differences between the Jastrow and 
UCOM corrected multipole contributions increase with 
increasing multipole. This pattern is reminiscent of the 
one shown in Fig. 3 of [l^l where no finite-size or higher- 
order-term corrections were taken into account. All this 
bears evidence of the fact that the rudimentary Jastrow 
method overestimates the effect of nuclear short-range 
correlations in Oz^/3/?-decay calculations. 

In Table [J we display the effects of the various cor- 
rections to the matrix elements M^^'^'> of ^^Ge and ®^Se. 
There we show in the second column the bare matrix 
elements (b.m.e.), then in the third column we show 
the b.m.e. corrected for the higher-order terms in the 
nucleonic current (A). In the fourth column we have 
added the nucleon finite-size effect (B) to the previous 
matrix elements (b.m.e.+A), and finally, in the last two 
columns, we have added to the previous matrix elements 
(b.m.e.+A+B) either the Jastrow (C) or UCOM (D) 
short-range corrections. 



TABLE II: Matrix element {M^^'^'^l of Eq. © 
tained in the present calculation and by Simkovic et al. |36l |. 
Shown are the results without and with the short-range cor- 
relations (s.r.c.) for gA = 1.254. 



without s.r.c. 
Present [3 6] 



Jastrow 



with s.r.c. 
UCOM 



6.36 



5.16 



4.72 



1.08 



TABLE III: Beta decay log ft values for transitions from the 
2^ states of ^®As and *^Br to one- and two-phonon states in 
the indicated final nuclei. 



= Se 



^Kr 



Jf 


exp. 


th. 


exp. 


th. 




9.7 


9.0 


8.9 


9.3 


2^ 


8.1 


7.7 


7.9 


7.7 




10.3 


9.2 


> 9.6 


9.4 




8.2 


8.7 


8.0 


9.0 




11.1 


10.9 


? 


11.1 



TableUis very interesting in the sense that there we can 
access the magnitudes of the various corrections to the 
bare matrix element. The magnitudes of the corrections 
coming from the finite nucleon size and the higher-order 
terms of the nucleonic current together amount to 25% - 
30%. In fact, the magnitude of our bare matrix element 
is rougly the one reported in j2ll,[33|- Even after correct- 
ing by the higher-order terms and the nucleon finite-size 
effect the matrix elements agree roughly, as shown in Ta- 
ble ini From the table we also see that our computed 
Jastrow corrections are much less than the ones obtained 
in jU, llHl- The reason for this is not clear. The net 
effect is that our Jastrow, and especially the UCOM cor- 
rected matrix elements are much larger than the ones of 
[m, [IBl- On the other hand, our shell-model computed 
Jastrow correlation corrections for "^^Ca js^] agree with 
the ones of [4l| . 

To check the consistency of our calculations we also 
computed the single /?~ decay rates from the lowest 2~ 
states of the intermediate nuclei ^^As and ^^Br to the 
lowest 2+ collective state, 2^ , and its Oj, 2^ and 4]*" two- 
phonon excitations (see, e.g., [3^) in ''^Se and ^^Kr. The 
wave function of the 2^ state was calculated by the use of 
the quasiparticle random-phase approximation (QRPA) 
[l^ l and its energy was fixed to the experimental one [3^ . 
Beta decays to the mentioned final states were computed 
by the method of the multiple commutator model (MCM) 
of H. 

The 2~ wave function was calculated by using the cen- 
tral value of gpp in the intervals (0]). In the case of a 
2~ initial state this choice works fine since the calculated 
beta-decay rates depend only weakly on gpp within the 
range relevant for the 2vl3l3 and Oi^/3/3 decays. We com- 
pare the computed log/i values with the available data 
in Table IIIII From this table it is seen that the com- 
puted numbers nicely reproduce the trends of the mea- 
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TABLE IV: Nuclear Oi//3/3 matrix elements of Eq. ^ for the 
decays of ^^Ge and *^Se. The UCOM and other corrections 



are included. 


The used Qpp 


values 


are also indicated. 










«2Se 




9pp 


1.02 


1.06 


0.96 


1.0 




1.923 


1.803 


1.304 


1.214 


M^V 


-4.632 


-4.208 


-3.293 


-2.950 




-6.555 


-5.355 


-4.597 


-3.722 



sured ones, although the assumpion that the states Oj, 
2^ and 4^*" are pure two-phonon excitations is an ideaHzed 
one. It has to be noted that there is no experimental data 
on beta decay of the lowest l"*" state in ''^As and ®^Br. 
For some double-beta-decaying systems the 1"^ data ex- 
ists and matching of beta and double beta decay could 
be more problematic due to the stronger gpp dependence 
of the beta-decay rates from a 1+ intermediate state [25| . 

Our final results for the OvP/S nuclear matrix elements 
have been collected in Table [iVl These matrix elements 
were calculated with the UCOM short-range corrections 
by taking into account also the finite-size of the nucleons 
and the higher-order terms in the nucleonic weak cur- 
rent. The two different values of the matrix elements 
correspond to the gpp and gA parameter combinations 
indicated earlier in the text. As can be seen, the values 
of the final matrix elements vary between 

5.355 < IM^O"^)] < 6.555 for ^^Ge , 
3.722 < |m(°'')| < 4.597 for ®^Se . 

We compare these matrix elements with other recent 
calculations in Table |Vl There the values 1.0 < gA < 
1.254 are used for the axial-vector coupling constant, ex- 
cept that Simkovic et al. [s^ use = 1.254. The results 
of Civitarese et al. [1^ are based on the formalism in- 
troduced in [i^l where the finite-size of the nucleon and 
the nucleonic weak current were obtained from a rela- 
tivistic quark-confinement model. In ji^l the generated 
nucleonic weak current is incomplete as compared to the 
present formalism, adopted from |36i]. Also the short- 
range correlations were not taken into account. In this 
sense the last column of Table [V] should be compared to 
the fourth column "+ A+B" of Table HI As already said, 
the differences between the two results can be explained 
by the different treatment of the weak- interaction current 
and the nucleon form factor. The above stated, we con- 
clude that our present results are more complete that the 
ones of [2^ and thus should be more reliable. It is worth 
mentioning that in [s^ the quoted matrix elements were 
calculated for the "default" value gpp = 1.00 without tak- 
ing into account the higher-order terms in the nucleonic 
current and the nucleon finite-size effect. 

In the matrix element calculations there may be other 
uncertainties than the ones induced by the uncertainty 
in the value of gpp. Such uncertainties could arise 
from sources such as deformation, the mean-field single- 



TABLE V: Values of the matrix element \M^°'''^\ of Eq. © 
obtained in several recent calculations. 



Nucleus 


Present 


[21] [36] 


[26] 


''^Ge 


5.36-6.56 


2.26 - 2.74 2.80 


4.03 - 5.92 


«2Se 


3.72-4.60 


1.86 - 2.45 2.64 


2.82 - 4.14 



particle energies, and the adopted two-body interaction. 
In [21I it was shown that the effect of the adopted two- 
body interaction is very small as long as the interaction 
is microscopic. Our adopted Bonn-A interaction is of 
this type and included in the survey of [2l|. In [2l| it 
was furthermore demonstrated that the size of the single- 
particle space does not produce sizable effects as long 
as the value of gpp is determined from the 2vj3(3 data, 
as is done in the present calculations. By the same ar- 
gument, only small effects are expected from different 
parametrizations of the Woods-Saxon mean-field poten- 
tial and the resulting slightly different single-particle en- 
ergies. These two sources of uncertainty produce effects 
that can be expected to be smaller than the one coming 
from the short-range correlations. 

The role of deformation is the most uncertain one. The 
presently discussed nuclei of the two double-beta-decay 
chains are pf-shell nuclei and most likely they have no 
or very small static deformation. Instead, they are most 
likely soft anharmonic vibrators. The deformation allows 
of a new suppression mechanism of 2v(3j3 decay, namely 
through the overlap factor used to take into account the 
non-orthogonality of the intermediate states generated by 
using the initial and final ground states as starting points 
in two separate pnQRPA calculations. This suppression 
mechanism is enhanced when the deformations of the ini- 
tial and final nuclei of double beta decay are different 
[i^ j. However, in [i^l it was deduced experimentally that 
^^Ge and ^^Se exhibit quantitatively very similar neutron 
pair correlations. This would indicate similarity of their 
ground states and no suppression would arise through dif- 
ferent ground-state deformations. For ®^Se and *^Kr this 
question is still open. In any case, the role of deformation 
in Qvj3[3 decay is still largely unexplored and no definitive 
conclusion about the importance of deformation can be 
drawn for the present. 

Our final matrix elements can be converted to Qv(3l3 
half-lives by choosing a value for the effective neutrino 
mass in ([2]) . Expressing the effective mass in units of eV 
and using the phase-space integrals tabulated in we 
obtain for the predicted half-lives 

4/2^ = (0.96 - 1.44) X 1024yr/((TO,)[eV])2 for ^^Ge , 

4/2' = (4.53 - 6.90) X 1023yr/((m,)[eV])2 for ^^Se . 

(6) 

In summary, we have calculated the nuclear matrix ele- 
ments for the Qvfif] decays of ''^Ge and ^^Se by using the 
proton-neutron quasiparticle random-phase approxima- 
tion with a realistic two-body interaction and a realistic 
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single-particle space. The numerical calculations were 
done by including the higher-order terms of the nucle- 
onic weak currents, the nucleon's finite-size corrections 
and the nucleon-nucleon short-range correlation effects. 
The short-range correlations have been calculated by us- 
ing the unitary correlation operator formalism that is 
superior to the traditionally adopted rudimentary .Jas- 
trow procedure. The UCOM reduces the bare values 
of the computed matrix elements less than the Jastrow 
procedure, leading to larger matrix elements than the 



ones quoted in the recent literature. This reduces the 
predicted theoretical half-lives of ''^Ge and *^Se 

by a significant amount and thus directly infiuences the 
neutrino-mass sensitivities of the running and future dou- 
ble beta experiments. 
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